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Observation of excitonic resonances in the second harmonic spectrum of M0S2 
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We measure the second-harmonic (SH) response of M 0 S 2 in a broad fundamental photon energy 
range between 0.85 and 1.7 eV, and present a continuous SH spectrum capturing signatures of both 
the fundamental A and B excitons in addition to the intense C resonance. Moreover, we interpret our 
results in terms of the first exciton simulation of the SH properties of multi-layered M 0 S 2 samples, 
represented by a study of trilayer M 0 S 2 . The good agreement between theory and experiments 
allows us to establish a connection between the measured spectrum and the underlying electronic 
structure, in the process elucidating the non-linear excitation of electron-hole pairs. 


In recent years, two-dimensional transition metal 
dichalcogenides (TMDs), being semiconductors 1 4 I with 
band gaps in the 2-3 eV range 4 ®, have been applied as 
the basis for a host of novel two-dimensional optoelec¬ 
tronic devices® 3 . Accordingly, a great deal of attention 
has been devoted to the linear optical properties of these 
compounds 1 3 . Also, non-linear optical techniques have 
been demonstrated as particularly powerful probes for 
the microscopic structure of few-layered TMD c rystals, 
revealing e.g., their crystallographic orientation 7 23JELD 
and stacking angle 1 ^. Moreover, second-harmonic (SH) 
spectroscopy allows for insights into the electronic struc¬ 
ture of TMD flaked^, exposing edge-localized effect 
or valley-coherent excitation^. 

While experimentally determined linear response func¬ 
tions of various TMDs, domin ated b y electron-hole pairs 
bound by several hundred meVd^ are reproduced rea¬ 
sonably well at the Bethe-Salpeter equation (BSE) level 
of complexit}®®^, similar comparisons for non-linear 
cases are lacking. Although several proposed theoreti¬ 
cal mod els, ba sed on both independent-particle! 19 2 3 an( f 
excitonic 22 23 approaches, have been published, the SH 
spectrum of M 0 S 2 has been investigated experimentally 14 
only in a narrow fundamental photon energy range be¬ 
tween 1.2 and 1.7 eV, probing the so-called C resonance 
also known from linear optics. Theoretical models have 
so far been benchmarked by their ability to reproduce 
this single feature 22 —. However, the parabolic disper¬ 
sion, and the split valence bands near the K-points of 
the Brillouin zone, translates into bound electron-hole 
pairs in the BSE picture. These give rise to sharp peaks 
in the absorption spectrum due to the so-called A and B 
excitons 3 ^ at photon energies near 1.9 eV, and the ques¬ 
tion remains how they might affect the SH response. In 
particular, the relative intensities and lineshapes of such 
features are of interest. 

In the present letter, we report an experimental op¬ 
tical SH spectrum generated from many-layered M 0 S 2 , 
with fundamental photon energies varied in a broad range 
between 0.85 and 1.7 eV, capturing both the fundamen¬ 
tal A and B excitons in addition to the aforementioned 
C resonance. We perform our optical experiments on 
(i) many-layered M 0 S 2 flakes exfoliated onto fused sil¬ 
ica by the well-known scotch tape method^, and (ii) a 
bulk M 0 S 2 crystal. The first approach gives access to 
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Figure 1. (a) Schematic representation of N 2H-stacked M 0 S 2 
layers, each with a sheet SH susceptibility of alternating sign. 
An incoming fundamental intensity I u induces SH reflected 
R ^ and transmitted T ^ signals, (b) Microscope image of 
thick M 0 S 2 flakes exfoliated onto fused silica. The image is 
140x140 11 m. (c) Atomic force micrograph of a few selected 
M 0 S 2 flakes. The inset displays height (z) profiles along the 
lines indicated on the main figure. 


a relatively smooth surface suitable for optical experi¬ 
ments, with little diffuse scattering. However, here flake 
thickness and surface coverage are difficult parameters to 
control within the area of a mm-sized laser beam. The 
thickness of the bulk crystal, on the other hand, is im¬ 
material due to absorption. Unfortunately, the rather 
rough surface of weakly van der Waals bound bulk TMD 
crystals typically make systematic optical studies diffi¬ 
cult. Hence, due to large amounts of diffuse scattering, 
we only measure reliable results for maximized SH sig¬ 
nals at the C resonance. Thus, the flake samples give 
access to highly resolved SH spectra, whereas absolute 
values may be extracted from the bulk sample, which, in 
turn, yield very little spectral information. By combining 
the measured SH response from these two sample types, 
we are able to calibrate the intensity measured from the 
exfoliated flakes against the bulk crystal measurements, 
thereby allowing for estimation of absolute parameters. 

In Figs, [ljb) and (c), we include microscope images of 
the produced flakes. Atomic force microscopy and pro- 
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filer scans confirm most flakes to consist of 40 - 100 mono- 
layers (MLs). Indeed, no measurable photo-luminescence 
(PL) was observed, regardless of pump wavelength. This 
suggests very little ML coverage, since only ML domains 
posses the direct gap necessary for efficient PL 1 3 . 

At first glance, dipole selection rules prevent genera¬ 
tion of even-ordered non-linearities in M 0 S 2 due to in¬ 
version symmetry, which is only broken in finite sam¬ 
ples consisting of an odd number of 2H-stacked M 0 S 2 
layers 14 16 . On the other hand, thick flakes, such as those 
investigated here, are not expected to display properties 
markedly different from bulk samples. From a simplistic 
point of view, the layer-selective vanishing SH response 
may be understood from the fact that adjacent layers are 
mirror images of each other in the yz- plane, as depicted 
in Fig. [lja). Hence, in the limit of weak inter-layer cou¬ 
pling, adjacent layers contribute anti-symmetrically to 
the SH response function While this means that 

the net polarization of two adjacent M 0 S 2 layers can¬ 
cels in the dipole approximation, the spatial variation 
of the driving field makes this cancellation only approx¬ 
imate. Indeed, previous measurements on bulk M 0 S 2 
crystals confirm this, with non-vanishing (although very 
weak) SH signals reportecP^U. It is this effect that 
allows us to perform SH spectroscopy on many-layered 
M 0 S 2 samples. We here analyse the measured response 
in terms of an effective, spatially dependent SH suscep¬ 
tibility x^ 2 \ z )-> constructed from the effective sheet sus¬ 
ceptibilities ix^eff the individual MLs, viz 

N 

X (2 \z,lo) = -Xsl s (u)^2(-l) n 5{z - nd + d). (1) 

n= 1 

Here, d ~ 6.2 A is the centre-to-centre interlayer distance 
depicted in Fig. [lja). 

The SH experiments are performed using a Nd:YAG 
pumped optical parametric oscillator, yielding 5 ns pulses 
at a repetition rate of 10 Hz and average pulse ener¬ 
gies of approximately 1 mJ at the sample. Measure¬ 
ments on flake and bulk samples are recorded in trans¬ 
mission and reflection configuration, respectively, using a 
small optical angle of incidence approximating 4°. Fur¬ 
thermore, the spectral distribution of the emitted ra¬ 
diation was analysed with great care to rule out con¬ 
tributions from competing non-linear processes, such as 
two-photon PL. All results are reported relative to the 
SH intensity reflected from an a-quartz (Qrz) wedge, re¬ 
lated to the pump intensity I u by the expression I^ z = 
I^Qrzl 2 IXQr Z | 2 ^2- The bulk SH susceptibility of quartz 
has been measured tcP^ Xqi ~ 0.3 pm/V, while the SH 

reflection coefficient p^ vz has been derived beford^. 

We also note that the applied spot size is on the 
order of 1 mm 2 , hence, encompassing multiple crys¬ 
tal domains of random orientation. For this reason, 
we are not able to observe rotational anisotropy, which 
is routinely seen when investigating few-layered TMDs 
by SH microscopy 14 16 . Further, the various domains 
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Figure 2. SH signal relative to a quartz reference for thick 
M 0 S 2 flakes deposited on fused silica. The inset displays the 
emission spectra for a selection of pump photon energies. No 
two-photon PL is observed at photon energies corresponding 
to half the M 0 S 2 fundamental exciton energy. 

are mapped onto separate areas of the PMT detector, 
which means their respective contributions add incoher¬ 
ently, ruling out interference effects from adjacent crystal 
grains. 

In Fig. [2j we present the measured SH transmission in¬ 
tensity spectrum I ^ for a sample with M 0 S 2 hakes. We 
observe clear maxima in the SH spectrum at fundamental 
photon energies corresponding to half the energies of the 
linear absorption features denoted A, B, and C in Ref. [25l 
Similar maxima appear in our simulation 22 of the exci- 
tonic SH response of ML M 0 S 2 , where the corresponding 
peaks in the SH spectrum were dubbed A/2, B/2, and 
C/2, as also indicated on Fig. |2j Such ML simulations 
do not include effects of interlayer coupling, which are 
obviously present in the samples investigated here. To 
consider effects of multiple layers, we have extended our 
exciton modei 22 ^ to study SH generation in trilayer (TL) 
M 0 S 2 , which represents a minimal M 0 S 2 multilayer hake 
without inversion symmetry. We expand exciton states 
in a basis of singly excited Slater determinants, which 
are, in turn, constructed from tight-binding orbitals de¬ 
rived from the parameters of Ref. 29. For MLs, the 
low-energy spectrum is dominated by transitions from 
the spin-orbit split top valence bands to the degenerate 
conduction bands near the K-points of the Brillouin zone, 
giving rise the A and B excitons, as mentioned. In bi¬ 
layer M 0 S 2 , a very similar picture with two split, doubly 
degenerate top valence bands and a degenerate conduc¬ 
tion band minimum near the K-points is observed. Here, 
however, the splitting is not a spin-orbit effect, but due 
to crystal held splitting arising from interlayer coupling. 
In TL M 0 S 2 , we observe a combination of spin-orbit and 
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crystal field splitting, giving rise to two groups of triply 
near-degenerate valence bands, which preserve the low- 
energy A/B exciton structure. When adding more lay¬ 
ers, transitions in the energy range corresponding to the 
C peak become dominated by additional bands, which 
are not degenerate, resulting in a broadened C-feature. 
See the supplementary material 30 for more information. 
We apply a relatively dense 60 x 60 &-grid, including 
24 valence and 24 conduction bands. Hence, this prob¬ 
lem is computationally difficult, and certainly beyond the 
scope of direct diagonalization. However, by carefully 
block diagonalizing the BSE based on out-of-plane sym¬ 
metry, and by using our Lanczos/Haydock approach for 
calculation of SH response functions^, we are able to 
generate converged optical spectra. Furthermore, as de¬ 
scribed in Ref. EE linear spectral features are increasingly 
broadened with increasing photon energy due to electron- 
phonon scattering. Here, we approximate this effect by 
implementing a frequency-dependent phenomenological 
broadening function, which we take to increase linearly 
with SH photon energies larger than the fundamental ex¬ 
citon energ}®. Note, the simulated TL sheet suscep- 
tibility Xs J * s calculated in the dipole approximation, 
hence, it represents the sheet susceptibility averaged over 
the three-layer structure while including interlayer cou¬ 
pling. Hence, this quantity is fundamentally different 
from XsL, which is an effective ML property valid for 
electronically decoupled layers. 

In Fig. [3j we include the calculated SH response func¬ 
tion for TL M 0 S 2 , and compare it with the square root 
of the observed SH signal. The latter quantity is related 
to the SH response function via | 2 |x^ | 2 /^, 

where / is the surface coverage and rffi is a SH trans¬ 
mission coefficient for an N layer flake 30 . As already 
discussed, both / and N are unknown quantities for the 
flake samples. However, by assuming rffi to depend only 
weakly on frequency, the square root of the observed SH 
signal may be compared directly with the SH response 
function |. Indeed, we observe the important A/2, 
B/2 and C/2 features both in experiments and theory, 
with similar relative intensities. The theoretical obser¬ 
vation of a C/2 splitting absent in the experiments may 
be attributed to the simplistic implementation of scat¬ 
tering effects. In f act, the linear case is often subject to 
similar deficiencies 6 22 . In the inset of Fig. [3J we com¬ 
pare ML and TL results. Although similar trends are 
found, we clearly observe an effect of including multi¬ 
ple layers, which is more pronounced than for the linear 
case^. For the TLs, we observe a slight enhancement of 
the A/2 and B/2 features relative to the ML case. Fur¬ 
thermore, in the studies using SH spectroscopy Malard 
and co-worker^ have demonstrated how the C/2 peak 
in ML samples is slightly red-shifted and broadened in 
TL samples (we include the spectra of Malard et a/P 4 in 
Fig. ! for comparison). The theoretical results shown in 
the inset reproduce this trend. 

We note that the observed SH intensity, displayed in 



Figure 3. Modulus of the SH response function compared to 
the square root of the measured SH transmission spectrum. 
We include also the experimental Ixs^l spectra for MLs and 
TLs published by Malard and co-workers 14 . Note, we scale 
these spectra by a factor 5 to fit them on the same figure. In 
the inset we compare our new TL calculation with our previ¬ 
ously reported ML spectrum 22 , however, now using frequency 
dependent broadening. 


Fig. i is roughly five times larger than the quartz ref¬ 
erence signal, even though the sample was covered by 
roughly 10% M 0 S 2 hakes. Hence, this hints at the ex¬ 
traordinarily larg^ 4 16 SH response of few-layered M 0 S 2 . 
Unfortunately, we are not able to extract absolute values 
of the SH susceptibilities from measurements on hake 
samples due to uncertain surface coverages and hake 
thicknesses. However, our measurements on a ~ 5 mm 
thick bulk M 0 S 2 sample are not affected by thickness due 
to absorption. For fundamental photon energies near 1.4 
eV, i.e. at the C/2 resonance, we measure a SH intensity 
rehected the bulk sample approximately 60 times larger 
than the SH quartz reference signal. To extract an effec¬ 
tive sheet susceptibility from this value, we note that the 

intensity of the rehected SH signal /JL ma y related 

( 2 ) 

to the fundamental intensity and by the rela- 

tion = \p { n\ 2 \x ( sIk\ 2i I, where is a SH reflection 
coefficient valid for an N -layer structure. The rehection 
( 2 ) 

coefficient p y N , suitable for an inhomogeneous non-linear 
response given by Eq. [lj may be found using a relatively 
compact expression^, which in the bulk (N —)> 00 ) and 
ML (N = 1) limits yield 


P (2) 

roo 


( 2 ) 

Pi 


Aiko 

[1 +n 2u ][l +n u } 2, 

8iko 

[1 +n 2w ][l +n u } 2 ' 


( 2 ) 

(3) 


Here, ko = uj/c while the optical phase shift upon prop- 













4 


agating through a single ML is taken to be negligible 
compared to the shift across all N layers. Furthermore 
n u and are, respec tively, the refractive indices of fused 
silica and MoSj^SU. We note that the ML limit Eq. [5] 
is identical to the response derive d starting from a non¬ 
linear surface (5-polarization^E3 and has been applied 
before for the study of few-layered TMD^ 4 Also, the 
bulk to ML intensity ratio at hu: = 1.4 eV may be eval¬ 
uated to \poo /P] 1 2 3 4 5 * * 8 ) | 2 ~ 0.15% , agreeing favourably with 
previously reported results 16 . It is also noted that the 
bulk expression Eq. [2] is very similar to the correspond¬ 
ing result for SH reflection from an infinite medium with 
homogeneous non-linear responsd^, with the main dif¬ 
ference being the ko proportionality, reflecting the fact 
that koXs ifi ac ^ s as an e ff ec ti ye bulk susceptibility. 

Since we measure the ratio 1 ^//^ rz = 

lx^effP~| 2 /IXQrzPQrzl 2 > and since we have expres- 
sions for both p^ and pQ^ z in addition to XqI z , w e may 
estimate an effective ML sheet susceptibility defined 

previously in Eq. 0 to lx.S’effl ~ nm 2 /V. 

Based on our measurements, we find the effective sheet 
SH susceptibility \Xs \ to be somewhat larger than the 
~ 0.1 nm 2 /V reported for M 0 S 2 MLs on quartz substrate 
in Refs. M and m On the other hand, in Ref. □2D sus¬ 
ceptibilities as large as 60 nm 2 /V are reported. In addi¬ 
tion to the experimental values mentioned, on-resonance 
moduli of x^ derived from theory are typically in the 
~ nm 2 /V range 22 211 Hence, the absolute value of the 
SH response function of M 0 S 2 is still a topic of consid¬ 
erable debate. We note that results generated using ad¬ 
vanced microscopy techniques for resolution of individual 
TMD flakes with well-defined thickness, and by applying 
well-documented reference measurements for extraction 
of absolute values, such as Refs. M and HU should be 
regarded with a high degree of validity. 


The bulk crystal applied in our work is composed 
of multiple grains, with different orientations. On the 
boundaries of such grains, the odd/even layer Xs an ti~ 
symmetry, implying dipole cancellation for perfectly 
stacked samples, does not lead to complete cancellation 
due to misalignment. Hence, these domains might give 
rise to an additional SH signal, perhaps partly explaining 
the large response measured here. 

In conclusion, we have extended the experimentally 
available SH spectrum of multi-layered M 0 S 2 flakes to 
include the broad fundamental photon energy range be¬ 
tween 0.85 and 1.7 eV. Moreover, we have performed an 
exciton simulation of the nonlinear optical properties of 
TL M 0 S 2 , and find the key features of the simulated spec¬ 
trum to be in good agreement with the experimentally 
observed spectrum. Hence, in this work, we have iden¬ 
tified experimentally the signatures of the A, B and C 
excitons, in addition to the relative intensities of the re¬ 
sulting peaks. Hence, these results extend the spectral 
region in which experimental data are now available for 
comparison with theoretically derived SH spectra, allow¬ 
ing for future comparisons with more advanced models. 
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